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Abstract 

Based on temperature and precipitation data from 15 meteorological stations and data on livestock populations, we digitized 
spatial data related to climate aridity and livestock density in Xilingol Leauge, a typical steppe region of northern China, and 
analyzed the relationships among vegetation cover change (as measured by NDVI), aridity and grazing activity. The results 
show that: (1) the spatial pattern of the relationship between NDVI and the de Martonne aridity index (Lm) was positive in most 
parts of Xilingol League; (2) livestock density increased from northwest to southeast, with 53.2% of the study area having 
livestock density of 40-80 standard sheep units/km2; (3) there was a binomial regression relationship between NDVI and 
livestock density, with higher livestock density that was concentrated in areas where the NDVI value ranged from 0.46 to 0.62; 
(4) there was complex linear correlation among NDVI, Lm, and livestock density; (5) partial correlation analysis revealed a 
significant negative correlation between livestock density and NDVI, and a positive correlation between Lm and NDVI, and the 
influence of aridity on NDVI was considerably greater than that of livestock density. 
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Introduction 

The Xilingol Steppe is located in central Inner Mongolia and forms part of the eastern Asia grassland subregion of 
the Eurasian grassland region. This natural grassland supports several abundant native vegetation and grassland 
types, and most of it is protected and intact. The grasslands and climate of the Xilingol Steppe support several 
sensitive and fragile ecosystems and provide ecological stability for this region of China. Climate change and 
intense human activities are impacting the grassland ecosystems in this region. As shown by earlier studies, these 
fragile grasslands are susceptible to damage and degradation if they are managed unsustainably, and over the last 
several decades the grasslands have suffered from degradation. 

In recent decades, Chinese and foreign researchers have conducted extensive research related to grassland 
degradation, showing that 18%, 39%, 60% and 73% of the grassland area in Inner Mongolia had experienced 
significant degradation by the 1960s, 1980s, 1990s and the 2000s, respectively. The fact that the grasslands of Inner 
Mongolia have experienced degradation is indisputable. 

In light of this situation, we used data related to de Martonne's aridity index (Dm), climatic conditions, and 
livestock density to characterize the influence of human activity on grassland vegetation cover, as represented by 
the normalized difference vegetation index (NDVI). 

Human activities in pastoral regions can be represented by livestock numbers, and the spatial distribution of 
livestock can reflect the effects of grazing on grassland ecology. Past research on the impact of grazing intensity on 
the health of grassland ecosystems has focused on the community and ecosystem levels, primarily using 
experimental plots to study the influence of variations in grazing intensity on soil, vegetation, microclimate and 
environment. 
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However, few research studies based on remote sensing technology have been designed to analyze the effects of 
human activities on grassland degradation. We digitized spatial data on livestock using Kriging's interpolation 
method, and referring to the population density interpolation method, we used the locations of residential areas 
and administrative divisions of towns as weight factors. 

The objectives of this study are to use meteorological and NDVI data along with statistical data on livestock 
distribution to examine how climatic variables and grazing activities affect grassland vegetation cover. Our aim is 
to gain the understanding of natural and man-made mechanisms that drive change in this grassland ecosystem. 
The findings of the study will be applied to develop scientific strategies and programs for grassland restoration, 
management and use patterns. 

Data and Study Methods 

Meteorological and Statistical Data 

Climatic data was obtained from 15 meteorological observation stations in Xilingol League (Figure 1), and 
statistical data on livestock numbers was obtained from county statistical reports. 
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FIG. 1 VICINITY MAP OF WEATHER STATIONS IN THE STUDY AREA. 


NDVI Data 

A Normalized Difference Vegetation Index (NDVI) was used as a measure of vegetation conditions in the study 
area. Remote-sensing data from 1981 to 2002 was extracted from composites of 10-day Advanced Very High 
Resolution Radiometer (AVHRR) data with a spatial resolution of 1.1 km. Sixteen-day MODIS NDVI data from 
2002 to 2007 with a spatial resolution of 250 m was also used. The NDVI values are modified to range between -1 
and +1 to make the analysis more convenient. MODIS NDVI data were re-sampled to 1.1 km to match with the 
spatial resolution of AVHRR NDVI data. The maximum value component was used to obtain the monthly and 
yearly data. 

Aridity Index 

An aridity index is an indicator used to characterizing the degree of aridity in a region, which in this case is 
characterized by the ratio of moisture and heat. De Martonne (1926) put forward a simple computational method to 
calculate an aridity index based on precipitation and temperature, using equation (1). 
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I dM = P/(T + 10). (1) 

Where IdM is de Martonne's aridity index, P is accumulated precipitation (mm) and T is average temperature (°C). 
Livestock Density Data 

Data on livestock numbers were converted into a uniform standard, using the standard sheep unit. Based on the 
1:250,000 scale maps and other administrative maps, we calculated the annual spatial distribution of livestock 
using the Kriging interpolation method. We also considered the locations of residential areas and administrative 
division of towns. That is, for each area, we calculated the number of villages and towns in the area and we defined 
the area of weight coefficient as 0.5 for these area. 


Results and Analysis 


Correlation between NDVI and Aridity 

Correlation analysis using ERDAS software was applied to the images of the spatial distribution of IdM and NDVI, 
producing a spatial distribution of the correlation index (Figure 2). There was a positive relationship between 
NDVI and IdM in the majority of the Xilingol area, except for in the Hunshandake Sandy Land and in the agro- 
pastoral transition zone in the southern part of the study area. Areas with correlation coefficients higher than the 
critical value of 0.31 accounted for 59.8% of the total area (P < 0.1), indicating that a significant positive correlation 
exists between moisture conditions and NDVI. 
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FIG. 2 CORRELATION BETWEEN ANNUAL NDVI AND DROUGHT INDEX (Idm) IN XILINGOL FROM 1981 TO 2007. 


Regression between Idm and Vegetation Cover 

Both vegetation cover and IdM increased in bands running from west to east across Xilingol League. The 
relationship between IdM and NDVI was also compared across three time periods representing the 1980s, 1990s, and 
the initial years of the 21st century (Figure 3). IdM was classified into different classes using 2.5 as an interval. The 
classification was done separately for each image from the three time periods. In each time period, IdM was derived 
for each class by averaging the IdM for that class Mean NDVI for each time period was also calculated by averaging 
NDVI in the same aridity index (Dm) class. 
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FIG. 3 THE RELATIONSHIP BETWEEN NDVI AND THE DROUGHT INDEX (Idm). 

Analysis revealed an obvious positive correlation between mean IdM and NDVI for all three periods (Table 1). 

TABLE 1 REGRESSION EQUATIONS BETWEEN NDVI AND ARIDITY INDEX (Idm) IN THREE TIME PERIODS. 


^~^Eqiiations 

Periocfs^^— 

Equation 

R2 

1980s 

y=2.4754x+8.9548 

0.9997 

1990s 

y=2.4572x+6.6365 

0.9996 

the beginning of 21st century 

y=2.4667x+6.4019 

0.9994 


Spatial Change in Livestock Density 

Livestock density for the period 1981-2007 increased in bands from the northwest to southeast. The majority of 
Xilingol (53.2%) supported a livestock density of between 40-80 standard sheep units/km2. Livestock densities 
lower than 40 standard sheep units/km2 occurred in Erlianhaote, the northwestern part of Suniteyouqi and the 
northern part of Suniteyouqi. Livestock densities higher than 100 standard sheep units/km2 occurred in the 
southwestern part of Xiwuzhumuqinqi and the southeastern part of the research region in the agro-pastoral zone 
(Figure 4 (a)). 

Compared with the 1980s, strong increases were seen in the 1990s in areas with a livestock density between 60-140 
standard sheep units/km2, while the areas of livestock densities lower than 40, between 40-60, and higher than 140 
standard sheep units/km2 decreased by 4.7%, 15.3%, 4.4%, respectively. Overall, livestock density increased in the 
1990s (Figures 4 (b) and (c)). In the initial period of the 21st century, livestock density decreased sharply and 
livestock density was under 120 standard sheep units/km2 across the entire region. Compared with the 1980s and 
1990s, livestock density of which areas lower than 40 and between 40-60 standard sheep units/km2 increased by 
11.9% and 16.5%, respectively (Figure4(d)). 

Livestock density is affected by a variety of human activities and by natural factors, including the climate. 
According to the analysis of the composite climate indicator (Dm), climatic conditions in the 1980s and 1990s were 
better than in previous periods, although the initial period of the 21st century experience relatively poor climatic 
conditions, with more extreme drought caused by higher temperatures and low precipitation. Climate affected 
vegetation cover, and herders most likely changed their stocking levels to respond to changes in vegetation 
conditions. The macro policy of the Chinese government also affected stocking levels. Grassland and livestock 
contracting systems were implemented in the 1980s in China, with use rights to grassland reverting to the villagers, 
and the market influencing herders' production decisions. The grassland tenure policy along with the use rights 
and contract responsibility system were implemented completely in the 1990s. With privatization of livestock, the 
number of livestock grazed increased greatly. The tenure, use right and contract responsibility system were mostly 
implemented in the early years of the 21st century, coupled with the implementation of transfer grazing and 
grassland enclosure policies. With the implementation of these policies, the use of grassland gradually transformed 
from extensive management to intensive management, and livestock density tended to be more reasonable. 


15 


www.seipub.org/fiee 


Frontiers in Environmental Engineering (FIEE) Volume 5, 2016 



FIG. 4 SPATIAL DISTRIBUTION OF ANNUAL LIVESTOCK DENSITY IN XILINGOL: (A) FROM 1981 TO 2007; (B) FROM 1981 TO 1989; (C) 

FROM 1990 TO 1999; (D) FROM 2000 TO 2007. 


Correlation between Vegetation Cover and Livestock Density 

Based on the ten-day NDVI data from 1981-2007, we obtained the average image for every decade and classified 
each image into different classes using five as the interval for each class. Livestock density was derived for each 
class by averaging the NDVI that belonged to the same class. Mean livestock density for each time period was also 
calculated for each NDVI class. Regression analysis was applied to analyze the relationship between mean NDVI 
and mean livestock density in each time period (Figure 5). 



FIG. 5 RELATIONSHIP BETWEEN LIVESTOCK DENSITY AND NDVI IN THE STUDY AREA 
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Analysis revealed a binomial regression relationship between NDVI and livestock density in Xilingol League. With 
an increase in NDVI, livestock density first increased and then decreased. High livestock density was concentrated 
in areas with NDVI values ranging from 0.46 to 0.62. Table 2 provides the regression equations. 


TABLE 2 REGRESSION EQUATIONS BETWEEN NDVI AND LIVESTOCK DENSITY DURING THREE TIME PERIODS 


1 ^uation 

Equation 

R2 


1980s 

y=-0. 5543x2+1 1 ,294x+25.626 

0.7647 

1990s 

y=-0.5936x2+12.765x+4.6742 

0.8676 

the beginning of the 21st century 

y=-0.6996x2+14.597x+18.946 

0.9300 


Mechanism of Vegetation Cover Change 

Using the vector data delineating the boundaries of villages and towns and the annual data on aridity, livestock 
density, and NDVI, we obtained a binary linear regression equation using ArcGIS software, as shown in Equation. 
( 2 ). 

z = 113.019 + 3.596x- 0.131y. (2) 

Where z stands for NDVI, and x and y stand for IdM and livestock density, respectively. The NDVI values are 
modified to range between 0 and 255 to make the analysis more convenient in this equation. 

The equation has an R2 of 0.772, F = 113.622 and P = 0.000. This equation revealed a complex linear correlation 
existing among vegetation cover (NDVI), IdM, and livestock density in the study region. The relationship between 
NDVI and moisture conditions is positive, while the relationship between NDVI and livestock density is negative. 
This indicates that vegetation cover is influenced by climate and overgrazing in Xilingol. 

In order to determine which factor influenced NDVI the most, we did partial correlation analysis to measure 
multiple variables to determine the degree of linear correlation between two variables while using the other 
variables as controls. This provided a reasonable and reliable way to describe the linear connections between two 
variables. 

The results of partial correlation analysis revealed a significant positive correlation between livestock density and 
IdM (with a partial correlation coefficient of 0.811), which indicates that a more mesic environment can result in 
higher livestock stocking densities. A significant negative correlation was found between livestock density and 
NDVI (with a partial correlation coefficient of -0.336), while a positive correlation was found between the aridity 
index IdM and NDVI (with a partial correlation coefficient of 0.706). Moreover, the influence of aridity on NDVI was 
considerably greater than that of livestock density. 

Discussion and Conclusion 

At present, research on the effects of grazing on changes in grassland vegetation cover has mainly been done using 
experimental tests in the study region. However, the spatial distribution of grazing is also influenced by the 
distribution of human settlements, the extent of livestock grazing activities and the policies in place in different 
locations. The distribution of livestock grazing is not uniform across the landscape, so at a large scale it will be 
difficult to explain spatial differences in the correlation between grassland conditions and grazing activities if 
conclusions are drawn based only on data from small-scale experiments. We were able to use spatial data on 
livestock populations with Kriging interpolation techniques, considering the number of villages and towns in an 
area to interpolate the spatial distribution of population density. Our analysis revealed how climate and grazing 
are driving changes in grassland condition, taking NDVI as a proxy for vegetation cover, in Xilingol League. 

(1) In 1883, Dokuchaev proposed theories related to soil zone theory and the zonal pattern of the distribution 
vegetation, which have been developed into the most important fundamental theories in ecology. These theories 
have profoundly demonstrated the close relationship between the spatial distribution of vegetation and climate at 
both regional and global scales. The zonal distribution pattern of vegetation reflects the effect of climate on 
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vegetation development and distribution that is based on a long history of natural geographical phenomena. But 
within the time scale of a few years, with the current background of ongoing climate change, will there be a 
corresponding change in the relationship between growing conditions in this steppe region and climate change 
induced drought? The results of this study revealed a significant positive correlation existing between the mean IdM 
and NDVI in three time periods. The spatial pattern of the relationship between NDVI and IdM was positive in the 
majority of Xilingol League. That is, NDVI increases as moisture conditions improve, and decreases as moisture 
conditions decline. 

(2) At present, the identification of the dominant factor causing change in grassland condition is still being debated. 
Hong Fuceng, the Group Chairman of the China Agricultural Expert Consulting China Grassland Honorary 
Chairman of the Society, thought that global climate change, in particular climate warming caused by the 
greenhouse effect, has exacerbated the current trend of widespread deterioration of the grassland environment. 
Also, human activity is another important factor affecting grassland degradation. Professor Ren Jizhou of the 
Chinese Academy of Engineering, Lanzhou University College of Pastoral Agriculture Science and Technology has 
said: "The causes of grassland degradation include natural and human factors, although human factors are the 
main factors" (2008). These comments reflect two aspects of the problem as seen by two experts, highlighting 
climate change and human disturbance as the main causes of grassland degradation. The results of this study can 
provide some insight into the possible effects of climate change-induced drought and grazing pressure on 
grassland degradation. There was a binomial regression relationship between vegetation cover and livestock 
density in Xilingol League. The relationship between NDVI and moisture conditions was positive, while the 
relationship between NDVI and livestock density was negative, indicating that spatial variation in vegetation cover 
in Xilingol is influenced by both climate and overgrazing. The results of partial correlation analysis revealed a 
significant negative correlation between livestock density and NDVI and a positive correlation between IdM and 
NDVI. Aridity had a stronger influence on NDVI than livestock density. 
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